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5 Abstract. In a seed plant a pollen tube is a vessel that transports male 

e gamete cells to an ovule to achieve fertilization. It consists of one elongated 

7 cell, which exhibits growth oscillations, until it bursts completing its func- 

s tion. Up till now, the mechanism behind the periodic character of the growth 

9 has not been fully understood. It is shown that the mechanism of pressure 

10 - induced symmetry frustration occurring in the wall at the perimeter of 

11 cylindrical and approximately hemispherical parts of a growing pollen cell, 

12 together with the addition of cell wall material, suffices to release and sustain 

13 mechanical self-oscillations and cell extension in pollen tubes. At the tran- 

14 sition zone where symmetry frustration occurs and one cannot distinguish 

15 either of the involved symmetries, a kind of 'entangled state' appears where 

16 either single or both symmetry(ies) can be realized by the system. We antic- 

17 ipate that testifiable predictions made by the model (/ oc \/P) may deliver, 
is after calibration, a new tool to estimate turgor pressure P from oscillation 

19 frequency / of the periodically growing cell. Since the mechanical principles 

20 apply to all turgor regulated walled cells including those of plant, fungal and 

21 bacterial origin, the relevance of this work is not limited to the case of the 
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22 pollen tube. 

23 Keywords: cell wall; equilibrium equation; frustration potential; growth 

24 tensors; Nicotiana tobaccum; plant cytomechanics 

2 5 1 Introduction 

26 1.1 General outline 

27 The pollen tube has become a widely used cellular model system. In addition 

28 to being the fastest growing plant cell, it features periodic oscillations of the 

29 growth rate that have attracted numerous attempts to model the process^} 

30 While recent models have increasingly incorporated biological features such 

31 as ion transport and intracellular trafficking, a simple feature with potentially 

32 significant impact has been overlooked in past approaches: geometry. We 

33 modeled the strain rates in the cell wall caused by turgor pressure depending 

34 on the different symmetries present in the pollen tube and found that a 

35 crucial area on the cellular surface of the pollen tube is characterized by 

36 so termed symmetry frustration. This area represents the transition zone 

37 between hemisphere-shaped apex and cylindrical shank. From a biological 

38 point of view this zone is crucial since numerous molecular landmarks of 

39 polar growth are present on one side of this zone and absent from the other. 

40 The model predicts that the transition zone undergoes local peaks in 

41 strain rate opening intriguing avenues for research focusing on the polarity of 

42 the growth process. Furthermore, we propose that changes between different 

43 symmetry regimes might be the mechanical underpinning of periodic changes 

44 in growth rate and shape observed during oscillatory growth. We believe that 

45 our model make an important contribution to the field of plant cytomechanics 

46 in general and pollen tube growth in particular. 

1 The origin of this oscillation is still unclear, though all hypotheses agree in that, the 
cell wall mechanics are essential to the oscillation (Chavarria-Krauser and Yejie, 2011). 
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47 1.2 Preliminaries 



48 Pollen tubes are rapidly growing plant cells whose morphogenesis is largely 

49 determined by spatial gradients in the biochemical composition of the cell 

50 wall. Pollen tube growth is a critical process in the life cycle of higher plants 

51 (Winship et al., 2010). It has garnered a lot of attention and is at the center 

52 of considerable controversy (Kroeger and Geitmann, 2011a). The pollen tube 

53 is the carrier of the male gametes in flowering plants. A controversy swirls 

54 around the modes of extension leading to periodicity in growth and growth 

55 rate (Winship et al., 2010). While some authors claim that hydrodynamics 

56 is the central integrator of pollen tube growth (Zonia and Munnik, 2007, 

57 2009, 2011) leading to growth oscillations, the other couple the periodicity 

58 in growth dynamics to the changes in the wall material properties (Winship 

59 et al., 2010, 2011; Kroeger et al., 2011). 

eo Pollen tubes are tip growing cells, which means that cell wall unidirec- 

61 tional expansion is confined to the apex of the cell. They display extremely 

62 rapid growth that can be also reproduced under in vitro conditions. All 

63 growth activity - delivery of new cell wall material and cell wall deforma- 

64 tion - occurs at the tip of the cell (Geitmann and Cresti, 1998): the tube 

65 is capped by an approximately hemispherically shaped dome, the apex, to 
ee which all growth activity is confined (Fig. 1A). The deformation is driven 
67 by the turgor pressure, a hydrodynamic pressure inside the cell. Interest- 
es ingly, pollen tube evolution displays characteristic oscillations in growth and 

69 growth rate (Plyushch et al., 1995; Hepler et al., 2001; Feijo et al., 2001). 

70 The pollen tube (e.g. Lilium longiflorum, Nicotiana tobacum) growth oscil- 

71 lation depends on many underlying phenomena, amongst others the ion and 

72 mass fluxes, wall mechanical properties, system symmetry and turgor pres- 

73 sure. In isotonic conditions (Zonia and Munnik, 2007) the average growth 

74 cycle period T is about 50 s, while upon shifts to hypertonic or hypotonic 

75 conditions it is about 100 s and 25 s, respectively. The latter produce oscil- 

76 lations with typical frequencies (/ = 1/T): hypertonic - 0.01 Hz, isotonic - 

77 0.02 Hz and hypotonic - 0.04 Hz. The longitudinal and transversal oscillation 

78 power spectrum of an individual Nicotiana tabacum pollen tube (Haduch and 
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79 Pietruszka, 2012) is visualized in Fig. 2, and can not solely be described by a 

80 double exponent model (Pietruszka, 2012; e.g. Fig. 8A), suitable for normal 

81 cell evolution, but incomplete for periodical growth. 

82 The cell wall is one of the structural key players regulating plant cell 

83 growth since plant cell expansion depends on an interplay between intra- 

84 cellular pressures and the controlled yielding of the wall (Geitmann, 2009). 

85 The cell wall may be treated as a polymer with the property of viscoelasticity 

86 (colloquially "elasticity"), generally having notably low Young's modulus e 

87 and high yield strain, which is a normalized measure of deformation repre- 
ss senting the displacement between particles in the body relative to a reference 

89 length, compared with other materials. Cell wall polymers are amorphous 

90 polymers existing above their 'glass' transition temperature, so that consid- 

91 erable segmental motion is possible. At ambient temperatures, the cell wall 

92 is thus relatively soft (e ~ 1 MPa) and easily deformed. Cells can grow to 

93 some extent simply by stretching their walls as they take up water. However, 

94 continued cell expansion involves synthesis of new wall material. Synthesis of 

95 cellulose at the plasma membrane and pectin and hemicelluloses components 

96 with Golgi apparatus deposits layers on the inside of the existing cell wall. 

97 A mechanical prerequisite for the unidirectional growth of pollen tubes for 

98 the (scalar) hydrostatic pressure is a softer cell wall at the tip of the cell, 

99 and more rigid at the distal part (Geitmann and Steer, 2006; Fayant et al., 

100 2010). This gradient of mechanical properties is generated by the absence 

101 or scarcity of callose and cellulose at the tip (Aouar et al., 2009) as well as 

102 by the relatively high degree of esterification of the pectin polymers in this 

103 region. The gradient in cell wall composition from apical esterified to distal 

104 de-esterified is reported to be correlated with an increase in the degree of cell 

105 wall rigidity and a decrease of visco-elasticity (Parre and Geitmann, 2005). 
we Also microindentation studies show (Fig. 2 in Winship et al., 2010) that 
107 the pollen tube tip is less rigid and that the distal stiffness may be opposed 
los to apical softness. Needless to say, to sustain growth processes, a balance 
109 between the mechanical deformation of the viscoelastic cell wall and the ad- 
no dition of new cell wall material must be achieved (Kroeger and Geitmann, 
m 2011b). 
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112 Turgor pressure is the pressure of the cell sap against the wall in plant 

n3 cells. This is a force exerted outward on a plant cell wall by the water and 

114 solutes contained in the cell. As a result it gives the cell rigidity. An ex- 

115 cess turgor pressure or cell wall local weakening can result in the bursting 
lie of a cell. Both constituents, turgor pressure and the wall properties are de- 
ny cisive for the mechanical properties and dynamics of the developing plant 
us cell. The osmotically maintained (hydrodynamic) turgor pressure in living 
n9 plant cells and the mechanical properties of the cell wall itself are among 

120 the most fundamental physical factors that dictate both cell growth and cell 

121 morphogenesis in plants (Schopfer, 2006). In fact, the physical properties 

122 of the wall and the turgor pressure have pivotal functions since they repre- 

123 sent the 'downstream parameters' of all cellular signaling events (Chebli and 

124 Geitmann, 2007). For our future purposes we note that turgor pressure is 

125 high in pollen tubes: 0.1 - 0.4 MPa in lily (Benkert et al. 1997; Winship et 

126 al, 2010). 

127 The pollen tube geometry can be described by two different symmetries - 

128 a hemisphere shaped apex, and a cylindrical shank zone, that are connected 

129 by a transition zone between the two parts. In the present paper we propose a 

130 mechanism of symmetry frustration occuring in this transition zone between 

131 the two involved symmetries as a possible mechanism responsible for growth 

132 rate oscillations. In simple terms by symmetry frustration we mean that a 

133 small ring of cell wall (hereafter referred to as interface T, Fig. IB) is unable 

134 to 'decide' if it should behave as an elastic cylinder or an elastic sphere. 

135 Following this hypothesis, oscillations may arise because this mesoscopic ring 

136 behaves as if it 'jumps' periodically between the two mechanical states of 

137 different capacity of strain energy. 

us The application of growth tensors to developing plant organs has been 

139 known for a long time (Kutschera, 1989). Such mathematical description 

uo has been utilized to apical meristems where the proliferating cells produce 

ui tissue stresses, which in turn influence the structure of the developing organ 

U2 and, hence the principal directions of growth (Hejnowicz, 1984). A differ- 

143 ent situation is observed for elongating plant organs, such as coleoptyles or 

144 hypocotyles, or elongation zones in roots, where cell division rarely takes 
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145 place. Also, the directional evolution of a single cell is primarily observed for 

U6 the elongating pollen tube. As it may be expected, various stresses occur in 

147 the elongating cell because of different properties of cell walls exposed to the 

us turgor pressure maintained by the gradient in the water potential (Kutschera, 

149 2000). The distribution of the wall stresses as well as deformation of the par- 

150 ticular wall layers can be calculated by solving equilibrium equations of elas- 

151 ticity theory. Such equations may help to locate deformation of a cell wall, 

152 exposed to an internal turgor pressure. The equilibrium equation may be 

153 derived both for materials deformed elastically (deformation vanishes when 

154 force equals zero) or non-elastically (plastic deformation survives, even when 

155 the acting force is removed). In this article we concentrate initially on elastic 

156 properties, because the highly controversial and uncertain mechanism of the 

157 oscillatory growth of pollen tubes is our main concern^] (Zonia and Munnik, 

158 2011; Kroeger et al., 2011; Winship et al., 2011). Though, plastic properties 

159 are inherently present in the proposed model via the derived (anharmonic) 

160 'frustration potential', and the assumed cell wall building processes located 
lei in the sub-apical, annular region, presumably at (about) the T - interface 

162 (Zonia and Munnik, 2008; Geitmann and Ortega, 2009). 

163 Possible mechanisms have been proposed to account for the oscillation of 

164 pollen tube growth rate in quantitative terms (Bartnicki-Garcia et al., 2000; 

165 Feijo et al., 2001; Dumais et al., 2006). A model for calcium dependent 
lee oscillatory growth in pollen tubes has been put forward (Kroeger et al., 2008). 
167 More recently the finite elment technique was used (Fayant et al., 2010) to 
lea establish biomechanical model of polar growth in walled cells. Also, the 

169 chemically mediated mechanism of mechanical expansion of the pollen tube 

170 cell wall by which deposition causes turnover of cell wall cross-links thereby 

171 facilitating mechanical deformation was set forth (Rojas et al., 2011). The 

172 role of wall ageing in self-regulation in tip-growth was considered in (Eggen et 

2 Whereas the hydrodynamic model as it is used by Zonia et al. proposes gradual 
increase in turgor until a threshold when rupture of individual links between cell wall 
polymers occurs, Winship and coworkers (Winship et al., 2010) state that turgor is es- 
sentially stable, but an exocytosis-induced relaxation of the wall causes expansion. They 
postuale that variations in cell wall mechanical properties cause the oscillations and that 
variations in turgor (if there are any) are a passive consequence due to cell wall relaxation. 
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173 al., 2011), while a model of plasma membrane flow and cytosis regulation in 

174 growing pollen tubes was discussed in (Chavarria-Krauser and Yejie, 2011). 

175 The irreversible expansion of the cell wall during growth as the extension of 

176 an inhomogeneous viscous fluid shell under the action of turgor pressure, fed 

177 by a material source in the neighborhood of the growing tip was examined 

178 in (Campas and Mahadevan, 2009). However, none of the models produced 

179 oscillations on mechanical basigS 

lso In our approach, which does not contradicts previous achievements, we 

lsi explore the relationship between turgor pressure and nontrivial cell geometry 

182 by changing loss of stability picture (Wei and Linthilac, 2007) to encompass 

183 cell wall mechanical properties in cylindrical and spherical geometries, both 

184 present in rapidly extending pollen tube. We base our physical model on the 

185 parametrized description of a tip growing cell that allows the manipulation of 

186 cell size, cell geometries, cell wall thickness, and local mechanical properties. 

187 However, the mechanical load (contrary to op. cit.) is applied in the form of 
las hydrostatic (constant) pressure. 

189 An important feature of pollen tubes elongation is that growth rate oscil- 

190 lates and, additionally, many of the underlying processes also oscillate with 

191 the same period, but usually with different phase (e.g. Fig. lc in (Zonia 

192 and Munnik, 2011)). However, the role of the oscillating ion gradients and 

193 fluxes in the control of pollen tube growth (Hepler and Winship, 2010) is 

194 beyond the scope of this paper and we will not discuss it here. Nonetheless, 

195 the outlined scenario leaves space for the periodical ion and mass fluxes in 
we and out of the cell. We also share the fundamental view expressed in (Rojas 

197 et al., 2011; Proseus and Boyer, 2006) that the wall extension is primarily 

198 a biophysical (mechanical) process, although ultimately dependent on enzy- 

3 The model presented by Chavarria-Krauser and Yeije (2011), actually does not de- 
scribe the cell wall, and hence, does not predict oscillations. The authors simply assume 
that the growth oscillations are given to understand the phase angle differences between 
growth velocity and the regulating mechanisms. However, by definition, oscillation is the 
repetitive variation, typically in time, of some measure about a central value (often a point 
of equilibrium), and pollen tubes do not exhibit such form of oscillations. In fact, what 
we observe in pollen tubes is a periodical elongation (without shrinking phase), and can 
be a result of transitions between two or more different states, which is another definition 
of oscillation, which we adopt in this paper. 
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199 matic activity, and that under conditions where the enzymatic background 

200 can be subtracted the biophysical process still proceedes normally. 

201 Any new model should deliver testifiable and quantitative predictions 

202 that can be validated by experimental data. In case of pollen tubes, it is 

203 necessary to present predictions that go beyond stress values which are in- 

204 herently difficult to measure. The presented model satisfies this requirement 

205 offering, among others, an experimentally testifiable power law (lu oc y/~P) 

206 between the turgor pressure P and growth oscillation angular frequency u. 

207 2 Results and discussion 

208 In search for the cause of experimentally observed pollen tube growth oscilla- 

209 tions we link analytic stress/strain relations with the mechanical properties 

210 of a tip growing cell, located at the perimeter of hemispherical dome and 
2n cylindrical shank. This is based on the observation that cell wall assembly 

212 by exocytosis occurs mainly at an annular region around the pole of the 

213 cell (Geitmann and Dumais, 2009; Zonia and Munnik, 2009) and that the 

214 concomittant turgor driven deformation of the cell wall causes characteristic 

215 strain in the hemisphere shaped apex of the cell (Fayant et al, 2010; Rojas 

216 et al., 2011). The dynamical properties of such a complex growing system 

217 should be found self-consistently (meaning that the turgor pressure and the 

218 wall mechanical properties are conjugate magnitudes that usually form cou- 

219 pled equations, which have to be solved by iteration methods). Nevertheless, 

220 in first approximation the following heuristic solution can be proposed. 

221 Assuming an intrinsic turgor pressure P, and a much smaller external 

222 pressure, of yet unspecified origin p ext (it can be just atmospheric pressure) 

223 producing an effective pressure P = P + p cxt ) the equilibrium equation for 

224 the displacement vector u, which is the shortest distance from the initial to 

225 the final position of a moving point, takes the form (Landau and Lifshitz, 

226 1986): 

2 (1 — v) grad(div u) — (1 — 2u) curl(curl u) = (1) 

227 where v is the Poisson coefficient, which is the ratio, when a sample ob- 
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ject is stretched, of the contraction or transverse strain, to the extension 
or axial strain. From now on Young modulus e(z) and Poisson coeficient v 
are assumed as picewise constant functions, so they remain constant on the 
interface T. 

Note, that by acting divergence operator on both sides of Eq. (1) we re- 
ceive A div u — 0, i.e. div u denoting the volume change due to displacement 
field is a harmonic function satisfying Laplace's equation. 

Eq. (1) may be solved analytically, providing that a problem exhibits a 
high degree of symmetry. In particular it can be solved exactly for spherical 
and cylindrical symmetries. Both symmetries are present in the description of 
pollen tube self-similar elongation, since the distal part resembles cylindrical 
tube while the apex is a hemispherical dome (Fig. IB). Thus, the symmetries 
present in both subdomains should be utilized in the description of pollen 
tube shape and dynamical properties. 

By assuming cylindrical symmetry (for the displacement vector field u = 
(iirjU^jUz) = (u r ,0,0), which is obtained under the assumption that the 
total length of the cylinder part remains constant (the axial elongation of 
the more rigid shank upon application of a constant internal pressure P we 
assume as negligible), and hence we accept u z = in the Ansatz (u r , u^, u z ) = 
(u r , 0, 0) instead of (u r , u^, u z ) = (u r , 0, u z ), which would lead to unavoidable 
numerical solution of Eq. (1)), and representing field operators (grad, div 
and curl) in cylindrical (polar) coordinates, Eq. (1) can be reduced to a much 
simpler form: 



This differential equation can be solved for the displacement field u r to yield 
the displacement for the cylindrical symmetry 



where a and b are constants to be determined from the boundary conditions 
(Landau and Lifshitz, 1986). Also, by introducing spherical coordinates with 
the origin in the center of a sphere the displacement field u is a function of 
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256 the radius r: u — (u r ,ue,u^) = (u r ,0, 0). Therefore rot u — and Eq. (1) 

257 reads: grad div u — 0. Hence, for the spherical symmetry we receive for 

258 displacement 

< = ar + ^. ( 4 ) 

259 where the upper index in Eqs (3) and (4) has been substituted to differenti- 

260 ate solutions for cylindrical (c) and spherical (s) geometries. The geometry 

261 of the elongating pollen tube can be described with a cylinder of radius tt 

262 capped by a half prolate spheroid with short radius tt and a long radius tl 

263 (see e.g. Fig. la in (Fayant et al., 2010)). Thus, as a model for normally 

264 growing tube we propose a thin-walled hollow cylinder ended by a hemispher- 

265 ical shell immersed in an external pool of pressure p ext and filled with a cell 

266 sap with turgor pressure P (we equate both radii tt = tl, for simplicity). 

267 The inner radius of the cylinder and a sphere is r±, while the outer radius 

268 is T2 (Fig. IB). Another simplifying assumption of the model is that we 

269 deal with weak (elastic) interactions at the interface T, due to wall building 

270 processes occuring at this region, and hence the deflection field may slightly 

271 differ on both sides (from the mechanical point of view they may be treated 

272 as weakly coupled). It is consistent with the view that deposition chemically 

273 loosens the wall by breaking load-bearing cross-links while simultaneously 

274 creating new, load-free cross-links, thereby effecting a fail-safe scenario for 

275 mechanical expansion (Rojas et al., 2011). 

276 Based on the proposition that the mechanical cell wall properties at the 

277 growing tip must be different from those in the shank, it was suggested 

278 (Dumais et al., 2004; Kroeger et al., 2008) that an anisotropy in the cell 

279 wall elasticity is required to account for the transition between spherical 

280 and tubular shape at the tip of the cell. Also, it was found (Geitmann and 

281 Parre, 2004) that the rigidity of the tip of the pollen cell was an increasing 

282 function of the distance from the apex. Therefore, the elastic properties 

283 of the cylinder imitating the cell wall at the shank and the hemispherical 

284 tip are represented by two pairs of material constants: Young's modulus e, 

285 also known as the tensile modulus, which is a measure of the stiffness of an 

286 elastic material and is a quantity used to characterize materials, and Poisson 
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287 coefficient v. In further calculations, we assume different values for these 

288 coefficients for distal (thick and rigid) and apical (thin and elastic) walls 

289 of a pollen tube. Such assumptions, about different values of mechanical 

290 constants at the apical dome and cylindrical shank, are fully justified (see 

291 e.g. Fig. 4 in (Geitmann and Parre, 2004) , where a spatial distribution 

292 of the Young's modulus is presented). Because we consider relatively small 

293 elastic deformations, the stress and strain tensors are related by the Hooke's 

294 law of elasticity (which is an approximation that states that the extension of 

295 a spring is in direct proportion with the load applied to it) and makes the 

296 deformation reversible. For the radial part of the stress tensor cr^ we have: 

297 a rr = —p cx t at r = r2 and a rr = — P at r = r±. Since the off-diagonal elements 

298 vanish, we are left with the strain u rr = ^ = a — \, = — = a + X and 

299 u zz = 0, the interesting radial a rr element of the stress tensor readsj 



300 By assuming boundary conditions as above, parameters a and b can be cal- 

301 culated (Lewicka and Pietruszka, 2009). They both depend on material con- 

302 stants: Young's modulus e and Poisson coefficient u, cylinder geometry (r 1; 

303 r 2 - radii) and pressure values P and p ext . 



305 Figs 3-4. In Fig. 3A we observe (a) lowering of deformation u r with 

306 radius r of apical hemispherical shell (b) almost constant u r for a thin-walled 

307 cylinder. This fact produces a wall stress presented in Fig. 3B equal to 

308 about 60 MPa (for the model parameters), which seems big enough to cause 

309 cell wall instability at T - interfac^J The calculated wall stress of the order 

4 In fact, to receive smooth solutions on T - interface equations for different geometries 
should be connected by transmission (gluing) conditions equating the forces and deflections 
on each side: an = a'n' and u = u', where n denotes the exterior normal to the boundary. 
In first approximation we let both subdomains be weakly coupled (visco - plastic phase) 
while cyclic wall building processes take place at T, and strongly coupled mechanically 
(visco - elastic phase) when wall building processes expire. We consider only the radial 
part of the stress a rr , deflection u r and strain u rr = d r u r tensors on T, for simplicity. 
This Ansatz, however, does not qualitatively influence the results. 

5 Typical plant cell turgor pressures in the range of 0.3 to 1.0 MPa translate into between 
10 and 100 MPa in the walls (Wei and Linthilac, 2007). 
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Quantitative calculations steming from Eqs (3) - (5) are presented in 
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310 of tenths of MPa is enough to cause local wall instability (radial strain), 

311 subsequent axial relaxation, wall building and unilateral cell expansion at the 

312 interface between the hemispherical apex and cylindrical shank. Unavoidable 

313 repetition of this process, owing to a constant effective pressure and a positive 

314 feedback mechanism necessary to drive oscillations (to overcome damping 

315 due to viscosity), may generate observable oscillations which continue until 

316 the wall building processes expire. Depending on the wall thickness, the 

317 calculated stress is equal to the distance between the curves (a) and (b) 

318 (inset). The tensile stress difference calculated at the apex and the distal 

319 part of the pollen tube cell wall (cr rr (apex) — a rr (distal)) is shown in Fig. 3C. 

320 Here it is parametrized by the turgor pressure P acting on the cell wall. The 

321 tensile stress in the wall clearly rises, as we increase the turgor pressure - 

322 this would, in turn, cause an increased growth oscillation frequency, as it is 

323 observed experimentally in the transition from hypertonic, through isotonic 

324 to hypotonic conditions (Zonia and Munnik, 2007). 

325 The influence of turgor on the oscillation period, as predicted by the 

326 model described in (Kroeger and Geitmann, 2011a), is also in agreement 

327 with our results presented in Fig. 4A and their consequences: the higher the 

328 pressure, the higher the oscillation frequency. Tensile stress difference at the 

329 transition zone, parametrized by the wall thickness, is presented in Fig. 4B. 

330 From the expression for tensile stress difference due to the symmetry 

331 change (note, curvature discontinuity and stress singularity occur at the tran- 

332 sition zone) a c rT — <y s rr = <j„. (cylinder) — a rr (sphere) which equals 

3 1 9 9 ^ 9 1 

r\ r 2 T\ r 2 T\ 7i 

f 3 (f2^ — fl^) ^ (j*2^ — ^"lj ^2^ — '"l^ ^2^ — fl^ 

333 and from the opposite formula: a s rr — a c rr , it is shown that calculations per- 

334 formed for an infinitesimally narrow ring T, where both geometries (cylin- 

335 drical and spherical) intervene, a symmetry frustration - leading to oscilla- 

336 tions of the radial part of the stress tensor - may take place. 'Frustration' 

337 originates from the fact that none of the locally involved symmetries is dis- 

338 tinguished. On the other hand, the calculated from Eq. (6) (by evaluating 



g.S p 
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339 E = f a Tr dr) strain energy density reads: 

E ^ ri2r 2 2 ( 2r3 + 2r ( r ! 2 + rir2 + r 2 2 ) _ rir2 ( ri + r2 )l (7\ 

2r 2 (r x 4 + ri 3 r 2 — r^ 3 — r 2 4 ) 

340 The quasi-discrete energy levels E_ and E + (possesing, however, a small 

341 dispersion 5E = 0.0003 [energy u.]) presented in Fig. 5 are non-degenerate 

342 due to the existence of a constant turgor pressure P which leads to the 

343 observed splitting (see also upper inset). Still, since both levels originate 

344 from the symmetry change at the critical limit considered in this work, they 

345 can be attributed to the oscillations in the pollen tube growth functions. 

346 Thus, the resonating frequency of growth (growth rate) corresponds to the 

347 energy difference E + — E_ = 2E (since E + = —E-), which in turn is directly 

348 proportional to the turgor pressure P. Consequently, we may equate the 

349 transition energy 2E between the resonating levels (Fig. 5) with pollen tube 

350 oscillation frequency observed in experiments (Eq. (7) implies that if P = 

351 then the system exhibits no oscillations, which is exactly the case, see also 

352 the plot of the potential energy U(r) at r = ro in Fig. 6). 

353 Notwithstanding, we note that the considered effect is exclusively con- 

354 nected with geometrically induced stress in the wall which may be linked 

355 with symmetry frustration^ and one can express it in measurable units [Pa 

356 m]. Indeed, calculating definite integral over the function expressed by Eq. 

357 (6): f r 2 [...]dr (with ri — 5 and r 2 = 5.25 jum) we receive the strain energy 

358 density: E± = ±0.256 [MPa /zm]. Therefore, the difference of strain energy 

359 density between the two levels is about 0.5 MPa for micrometer length scales 

360 typical for the width (which is about 250 nano-meters) of the pollen tube cell 

361 wall. Such energy density may lead to oscillations which are observable not 

362 only in growth rates. From our model it can also be deduced that the apical 

363 geometry oscillates (due to deformation u r located initially at T, compare 

364 Figs 3B and 4B in Pietruszka et al., 2012) to produce the so called pearled 

365 morphology (Rojas et al., 2011). Relicts (residues) of such deformations at 



6 A spatially degenerate ground state will undergo a geometrical distortion (the alter- 
ation of the original shape) that removes that degeneracy, because the distortion lowers 
the overall energy of the whole complex. 
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366 T are traced in Fig. lc and 6a, c (ibid.) as crests smeared out on a distance 

367 A, in agreement with our modeQ However, what else draws our attention 

368 is that there is no sign of deflecions at a distance shorter than the tube ra- 

369 dius R. The latter observation further supports the main idea presented in 

370 this paper of specific role of the T - interface in initiating oscillations. By 

371 assuming, after Rojas et al. (2011), the value of the linear t> avg = 0.2 [/xm/s] 

372 of the elongating cell and taking the average oscillation period T = 50 s 

373 from Fig. 7 we receive the observed wavelength of about A = 10 [//m], which 

374 is a doubled value of the radius R, as it should be expected assuming the 

375 correctness of our approach. In this picture, the distance A comprise the 

376 local deflection/wall stress/stress relaxation/recovery through wall building 

377 processes for every period T. 

378 Transitions between the states of different symmetry are shown in Fig. 5. 

379 The system is pumped with energy to jump over the analytic discontinuity 

380 (energy gap) between the (hemi-) spherical and cylindrical geometry. As a 

381 positive feedback mechanism necessary to drive the oscillations (to prevent 

382 damping due to viscosity) the energy is absorbed (ATP-pumping (Rounds 

383 et al., 2011)) in the transition zone above the state E_ producing the exited 

384 state E + . Then the system returns (by spontaneous symmetry breaking, to 

385 reduce the energy of the overall system) to the lower symmetry (cylindrical) 

386 state E- stimulating axial expansion. Both transitions (up and down) close 

387 one growth cycle with the prediced transition's rate w cc y, where T is the 

388 period. The whole process is repeated at the expense of pressure P and 

389 ATP-energy needed for wall synthesis (exhibiting growth oscillations), and 

390 eventually expires or reaches critical instability (the cell bursts). 

391 In a mechanical anharmonic oscillator, the relationship between force and 

392 displacement is not linear but depends upon the amplitude of the displace- 

393 ment. The nonlinearity arises from the fact that the spring (here: cell wall) 

394 is not capable of exerting a restoring force that is proportional to its dis- 

395 placement because of stretching in the material comprising the wall. As a 

7 Such geometrical oscillations of the wave-length A will be obtained when frustration 
occurs, and the cylindrical and spherical symmetries will be present on T - contour inter- 
changeably; compare Movie SI and S2 in (Pietruszka et al., 2012). 
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396 result of the nonlinearity, the vibration frequency and amplitude can change, 

397 depending upon the system elements displacement upon pressure P. An ap- 

398 proximate derivation performed in Appendix 1 delivers the analytic form of 

399 the (dual) 'frustration potential', which is a sum of attractive and repulsive 
wo forces, possibly responsible for experimentally observed growth rate oscilla- 

401 tions in pollen tubes. The latter, which is given by Eq. (12) and visualized 

402 in Fig. 6, we describe here shortly: Pollen tube oscillations are trapped at 

403 the potential well about the equilibrium point r for the corresponding sym- 

404 metric (harmonic) potential. Oscillation frequencies and amplitudes of the 

405 anharmonic potential depend upon the turgor pressure values, see Eqs (6) 

406 and (7), as it is observed (e.g. Fig. 4 in Kroeger et al., 2011; Kroeger and 

407 Geitmann, 2011a)). Wall expansion is allowed by molecular separation (r - 

408 values) exceeding those of harmonic potential. Dissociation energy at zero 

409 potential level corresponds to system instability (burst at T). The constant 

410 turgor pressure, as taken from Fig. 5, induces the value of U(r), hence the 

411 frequency of the oscillation and its amplitude. The low-lying (trapped) values 

412 deliver high frequencies and small amplitudes, while the higher-lying poten- 

413 tial values - low frequencies and larger amplitudes of oscillations. Above the 

414 critical threshold (corresponding to 'zero energy' at the vertical scale) a bond 

415 breaking occurs and the pollen tube burst at the transition zone, or deliver 

416 male gametes completing its function. The lower plot represents only one 

417 branch (of the prevailing cylindrical symmetry) of the full frustration poten- 

418 tial; the second branch (above) is in 'dual' subspace, and the oscillations take 

419 place between both branches. For the negative values of the cylindrical sym- 

420 metry, as shown in the plot, the mechanism of symmetry breaking favorizes 

421 this 'lower order' (cylindrical) symmetry for cell extension. 

422 The presented in Fig. 6 frustration potential is a more convenient model 

423 for vibrational structure of wall constituing molecules than harmonic oscilla- 

424 tor potential, because it explicitly includes the effects of bond breaking and 

425 accounts for anharmonicity of real bonds in the extending cell wall. It is 

426 also responsible for the inherent instability at the T - interface of a grow- 

427 ing tube (and - in consequnce - polymer building process), which can be 

428 experimentally supported by the fact that the pollen tubes always rupture 
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429 at the transition zone where the radial part of the strain tensor is consid- 

430 erable (Pietruszka et al., 2012, Movie S3). The form of the potential also 

431 contributes to the long debate among plant physiologists about the elas- 

432 tic/inelastic extension of plant cell wall in simple terms: any departure from 

433 parabola centered around r will lead to plastic extension, corresponding to 

434 elongation growth (Fig. 6). In addition, the infinite potential barrier at low 

435 distance r values prevents the growing cell wall from shrinking at a given 

436 pressure P. 

437 In order to calculate the value of the resonance angular frequency u> we 

438 momentarily accept the approximate (classical) relation: Eocw 2 . Assuming 

439 P = 0.3 MPa, taking the approximate A constant from the fit (see Fig. 

440 7) we receive / = 0.09 Hz, a value which belongs to the observed frequency 

441 spectrum in pollen tube growth functions (0.01 Hz - ~ 0.20 Hz, (Haduch and 

442 Pietruszka, 2012) for tobacco, Fig. 2; and (McKenna et al, 2009) for lily). 

443 As an aside, we stress that the calculated from Eq. (7) resonance frequency 

444 satisfies a power law u oc y/P (see Appendix 2 for detailed derivation). The 

445 application of this important relation to the experimental data is presented in 

446 Fig. 7. It is easy to notice, that this relation (u = 2ttA\/P, or equivalently 

447 / = Ay/P, where A is a constant - connected with the wall mechanical 

448 properties - to be determined from experiment, and [P] = MPa), if inverted, 

449 can serve (after calibration) to estimate difficult to measure turgor pressure 

450 P values from easy to measure oscillation frequencies (or periods T). 

451 Furthermore, it is clear that the material properties of the cell wall in the 

452 apical region should not be homogeneous, and therefore a proper mechanical 

453 description of growth must involve a gradient in material properties from the 

454 apical to the distal region (Fayant et al., 2010; Eggen et al., 2011). It has 

455 been shown (Eggen et al., 2011) that the calculated "expansion propensity" 

456 as a function of the distance from the apex measured in units of the tube 

457 radius R (notation, as in (Eggen et al., 2011)) shrinks to an area near the 

458 apex. Closer examination reveals that the inflection point is located at about 

459 1 pollen tube cylinder radius R, a place where we perform our calculations. 

460 The latter statement means that the slope is the greatest at z ~ R (in axial 

461 direction). This, and the fact that the 'dilution' sector is shown (Fig. 2 in 
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462 Eggen et al., 2011)) exactly at the limit of the two considered axisymmetric 

463 zones, is consistent with the view of intense changes of the wall mechanical 

464 properties at the limit of the distal and apical part. The 'dilution' effect is 

465 caused in our model by a rapid surface expansion due to displacement u r of 

466 the T - interface (see also (Parre and Geitmann, 2005): Fig. 1 - a local dip 

467 in the wall stiffness appears at about 10 /im from the apex, a place where 

468 our calculations are performed; and Fig. 7 (2) showing the possible local- 

469 ization of the transition zone from the cylindrical to spherical symmetry). 

470 Corresponding radial strain may trigger exocytosis that results in delivery 

471 of new cell wall material which rejuvenates this area. However, we should 

472 note, what we observe is not only according to the mechanical properties 

473 gradient, but mainly due to the changing symmetry at this place and the 

474 analytical consequences (curvature discontinuity) of this fact (compare e.g. 

475 Eqs 3 and 4). Likewise, observations, together with the analysis of Floures- 

476 cence Recovery After Photobleaching (FRAP) experiments (Geitmann and 

477 Dumais, 2009; Bove et al., 2008) indicate that exocytosis is likely to occur 

478 predominantly in the same annular region (cf. (Geitmann, 2010), Fig. 1) 

479 where wall expansion rates are greatest. It is concluded in the same work, 

480 that tip growth in plant cells does not seem to happen exactly at the tip. 

481 Further supporting data is provided in (Zonia and Munnik, 2008), where the 

482 vesicle fusion with the plasma membrane was shown to approach to within 

483 2-5 /im distal to the apex. The observations that growth is mainly at an 

484 annular region around the pole are in accord with calculations presented in 

485 this article. Probing mechanical properties at the perimeter of cylindrical 

486 and spherical part may result in calculation of the local rate of exocytosis 

487 (the conventional explanation for this phenomenon, i.e. oscillating exocyto- 

488 sis rate (see e.g. Fig. 4A in (Cardenas et al., 2008), Fig. 2 in (McKenna 

489 et al., 2009), is widely accepted). This may be roughly estimated by taking 

490 the oscillation frequency from Fig. 7. The read off value: / ~ 0.02 — 0.03 

491 Hz is in accord with the main observed periodical mode (Zonia and Munnik, 

492 2007) in the longitudinal power spectrum of pollen tube oscillatory motion, 

493 and presumably may be equated with the rate of exocytosis and new cell 

494 wall assembly in Nicotiana tobaccum pollen tubes and temporal variations in 



17 



495 the secretion of cell wall precursors (Kroeger and Geitmann, 2011a). This 

496 periodical growth activity, in turn, could be related (among others) to the 

497 relaxation of the stress (invoked by the hydrostatic turgor pressure) in the T 

498 - interface, i.e. in close proximity to the advancing apex of the cell, probably 

499 at the vesicle delivery zone in the subapex, where the vesicles are released 

500 into the cytoplasm (Fig. 3 in (Kroeger et al., 2009)). In addition, even 

501 small stress/strain fluctuations at this narrow cylindrical ring, belonging to 

502 both adjacent zones, localized on its circumference could lead to macrosopi- 

503 cally observable change of orientation of this ring and consequently direction 

504 change of the elongating pollen tube. This is indeed the case - see e.g. Fig. 

505 IF in (Zonia and Munnik, 2008), and even more pronounced in Fig. 1 in 

506 (Calder et al., 1997). It seems that polar growth in pollen tubes is associated 

507 with spatially confined dynamic changes in cell wall mechanical properties 

508 (Zerzour et al., 2009; and Appendix 2). Furthermore, the time course of an 

509 experiment, showing very little change in turgor pressure during cell growth, 

510 was measured by pressure probe monitoring of growing Lilium longiflorum 

511 (Winship et al., 2010, Fig. la) and re-analysed in (Zonia and Munnik, 2011, 

512 Fig. 3b-c). Even though direct measurements fail to indicate large-scale tur- 

513 gor changes during growth, rapid small-scale pressure changes (jumps) are 

514 visible, which presumably^] may be caused by the change of orientation of the 

515 tilt angles of wall building cellulose microfibrils in subsequent growth cycles 

516 in the considered transition zone^J As stated (ibid.), the measured periodic- 

517 ity for pressure oscillations ranges from 12 s to 25 s, which is the same as the 
sis routinely reported for oscillatory dynamics in lily pollen tubes (McKenna et 

519 al., 2009; Zerzour et al., 2009) and approximately agrees with the calculated 

520 frequency (/ = 1/12 = 0.09 Hz) from our model (see next paragraph). Closer 

521 examination of Figs 3a-c in (Zonia and Munnik, 2011) reveals, in addition, 

522 a slight but steady diminishing of turgor (negative slope) which can be a 



8 There is, however, no evidence for cellulose microfibrils to be involved in small changes 
in turgor. These small jumps in turgor may be simply imprccisions in the measurement 
method. This method requires readjusting the meniscus in a pressure probe needle for 
each data point and the measurement is inherently associated with significant noise. 

9 The direction of maximal expansion rate is usually regulated by the direction of net 
alignment among cellulose microfibrils, which overcomes the prevailing stress anisotropy 
(Baskin, 2005). 
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523 consequence of cell volume expansion in every cycle. The latter observation 

524 may be associated with the passive role of the turgor pressure in pollen tube 

525 growth, at least at unchanging osmotic potentials. In conlusion, we agree 

526 with the view, that in the growth process the (main) energy supply is de- 

527 rived from turgor, while the growth rate and direction from the wall local 

528 properties (Winship et al., 2010). 

529 Expansive growth in a plant cell relies on the interplay between the inter- 

530 nal turgor and the forces in the cell wall opposing deformation. Which of the 

531 two parameters controls the dynamics of growth has been controversial in the 

532 case of pollen tube growth (Zerzour et al., 2009). The long-standing model 

533 of pollen tube growth considers that cyclic changes in cell wall properties 

534 initiate growth (Holdaway-Clarke and Hepler, 2003; Winship et al., 2010). 

535 A new model based on accumulating data from recent work indicates that 

536 oscillations in hydrodynamic flow and intracellular pressure initiate growth 

537 (Zonia 2010; Zonia and Munnik, 2011). Both models agree that once growth 

538 initiated, osmotic pressure drives cell elongation. We have shown that the 

539 rapid polar growth phases during oscillatory growth in pollen tubes may be 

540 preceded by a strain - induced softening of the cell wall at the brink of the 

541 apical and distal parts (r). We also showed that cellular turgor pressure does 

542 not need to undergo changes during these repeated growth phases to display 

543 periodicity in growth. However, turgor pressure still preserves an important 

544 role, together with the cell wall mechanical properties, in controlling the dy- 

545 namics of pollen tube growth by changing wall stress and hence oscillation 

546 frequencies in the different osmotic environments - the frequency of oscilla- 

547 tory pollen tube growth in our model can be altered by changing the osmotic 

548 potential value of the surrounding medium (Kroeger and Geitmann, 2011a). 

549 There are many observations of oscillations that could affect growth rates 

550 (such as wall material deposition and extracellular ion fluxes). However, even 

551 from purely mechanical calculations, performed at the boundary between the 

552 wall cylinder shell and hemispherical shell at the apex, the following picture 

553 for the sequence of events for the elongating pollen tube emerges. A given 

554 (constant) turgor pressure produces different strain at the apical and distal 

555 wall parts possessing various mechanical properties and different symmetries. 
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556 This phenomenon is especially important at the narrow interface between 

557 these neighbouring parts. Consequently, a localized elevated stress of the 

558 order of tenths MPa is invoked in the wall causing serious instability at the 

559 brink of both sections generated by symmetry frustration, and the cell wall 

560 relaxation (loosening) process takes place in order to reduce tensile stress. 

561 This initiates a wall building process (which is an implicit assumtion of the 

562 model) in meridional direction. However, the effective turgor pressure pro- 

563 duces strain at the similar location (in the co-moving - with the moving tip - 

564 reference frame) axially equidistant from the tip and the whole process/cycle 

565 repeats. This eventually results in time - periodicity in the growth dynamics 
see recognized in the literature as pollen tube oscillations. 

567 3 Final comment 

568 This article offers a nontrivial solution for a long - sought mechanism of 

569 pollen tubes growth oscillations, which is the subject of swirling controversy 

570 in the field. It is based on the phenomenon of symmetry frustration of the 

571 cell wall in the apical region. This simple physical mechanism results in 

572 anlytically determined asymmetric 'frustration potential', the appearance of 

573 the landscape of 'discrete' energy levels at (different) constant pressures and 

574 the aforementioned oscillatory growth comes from the transitions between 

575 them. Moreover, a scaling relation between the turgor pressure P and the 

576 angular frequency of the oscillations u is derived, which is represented by a 

577 power 1/2 - law (u oc \fP). Later on, this prediction is successfully verified 

578 against a real plant physiological experimental data. 
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Appendix 1 



589 (i) The local force equation of motion in the mechanics of continuous bodies 

590 (due to Cauchy) reads (Lubliner, 2006): 



djUij + Fi = 0, 



(8) 
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591 where i, j = 1, 2, 3, and Einstein's summation convention is used. By ignoring 

592 any possible shear stresses and assuming for the infinitesimally small V - 

593 interface ass = and a zz = 



(9) 



594 we receive a rr = a(r), which is the only matrix element that survives. From 

595 Eq. (8) we have 

J djaijdx 1 + J Fidx 1 = (10) 

596 so since f Fidx 1 = —U(x) is the potential energy then f dja^dx 1 = U(x). 

597 Hence in our case J d r a rr dr = a rr = U(r), and one gets: U(r) = a rr = a(r). 

598 (ii) By identyfying 

U{r) = a{r) (11) 



599 we may follow Eq. (6) to receive 



a (3 



Aa ± = ±^T^ + C 



(12) 



where a = P^P^-, (3 = P and C is a constant. Next, in order to find 
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601 r we write Aa'(r ) = to get 

3a 2/3 

-^ + ^ = ° (13) 

602 Hence 

T " = 20' (14) 

603 Since U = Act we can plot the 'symmetry frustrated' potential U(r), Fig. 6. 

604 (iii) Considering small oscillations around equilibrium r , we may introduce 

605 new coordinate p: r = r — p and by substituting it to Eq. (12) receive 

u{ - x) = ^ {i- x y~ri{i-xy (15) 

606 where x = p/tq. By expanding both fractions for small x, we finally get the 

607 form for the harmonic potential: 



60s Comparing the above equation with the classical oscillator potential (m — 1): 

609 

U(x) = l -u 2 x 2 + U (17) 

6io and using Eq. (12) we get u 2 oc P. 

en Hence, the pollen tube oscillation frequency at the limit of small oscilla- 

612 tions equals u oc y/P, where P is the turgor pressure (compare also with Fig. 

613 7, where the proportionality constant (A) is estimated from experiment). 

614 Appendix 2 

615 Pollen tube oscillation: local deflection/wall stress/relaxation/recovery 

616 

617 STAGES of one oscillation (mechanical view): 

618 1. Recovery through wall building at T - interface: visco-plastic process 

619 (elastic equations do not apply here, because of wall and mass pro- 
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620 duction; the system is merely plastic and both subdomains are weakly 

621 coupled from the mechanical point of view). 

622 2. Strain and deformation production on T (the equations apply). 

623 3. Wall stress production on both sides of T and the resulting strain en- 

624 ergy: visco-elastic process (equations apply). 

625 4. Elastic strain energy relaxation in one cycle to produce elongation of 

626 one wave length (compare Rojas et al, 2011, Fig. 1C, Fig 6A): the 

627 phenomenology applies / = A\/~P, finding confirmation in comparison 

628 with authors' (Haduch and Pietruszka, 2012) performed experiment, 

629 Fig. 7. The process repeats: next oscillation takes place (go to 1. to 

630 start another oscillation). 
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Figure 1: Nicotiana tobaccum pollen tube apical region. (A) Microscopic view (B) 
Schematic view: radii of curvature r\ and r 2 , turgor pressure P and the investigated 
partition into two distinct regions (a narrow transition zone) are indicated in the chart. 
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Figure 2: Density plots of the longitudinal and transversal power spectrum of Nicotiana 
tabacum pollen tube obtained from the raw experimental data, bias subtracted (Haduch 
and Pietruszka, 2012) by Fourier analysis, calculated by the power of the Nyquist criterion 
and Nyquist rate, for different osmotic environments (-1 corresponding to the hypotonic 
case, - isotonic case, and 2.5, 3 and 5 corresponding to 25, 30 and 50 mM NaCl in 
hypertonic conditions, respectively). A broad, narrowing valley at the centre of the lower 
plot is clearly visible - low frequencies seen for longitudinal modes are shifted outwards. 
Red colour indicate high intensity peaks. Interpolated by DAVE, developed at NIST 
(Azuah et al., 2009). 30 




Figure 3: (A) Displacement u r due to the effective turgor pressure P (P = 0.3 MPa, 
Pcxt = 0.05 MPa) acting on the cell wall as a function of the radial distance r from 
the pollen tube long axis. Different Young modulus e and Poisson coefficient v in both 
subsystems: (a) u r — u s r for a hemispherical apex: r\ = 5 /im, T2 = 5.25 /mi, v = 0.4, 
e = 0.2 [GPa] (b) u r = u% for a cylindrical distal part: t\ = 5 /im, r2 = 5.25 /im, v = 0.2, 
e = 1 [GPa]. (B) Tensile stress a rr due to the effective turgor pressure P acting on the 
cell wall at the position where (a) the cylinder (shank) joins (b) the hemisphere (apex) as 
a function of the radial distance r from the pollen tube axis. Radial stress discontinuity 
between the distal wall and apical wall is proportional to the distance between the curves 
(a) and (b). The calculated maximum wall stress reaches about 60 MPa for the simulation 
parameters. The strain energy leading to oscillations is proportional to the shaded area. 
(C) Tensile stress difference <7 rr (apex) — a rr (distal) (in MPa) at the apex and the distal 
part parametrized by the changing turgor pressure P acting on the cell wall: (a) P = 0.5 
MPa, (b) P = 0.4 MPa and (c) P = 0.3 MPa. 



31 



5.00 5.05 5.10 5.15 5.20 5.25 5.30 
r[jum] 




Figure 4: Tensile stress difference (A) <r^ r — a s rrl upper curves, and the opposite case 
a s rr — a° r , lower curves, calculated at the boundary zone between the approximately hemi- 
spherical apical and the cylindrical distal part of the growing pollen tube. Parametrisation 
by the turgor pressure P acting on the cell wall (upper plots): P — 0.1 (green), P = 0.2 
(blue) and P = 0.3 MPa (violet). Remaining parameters for the respective wall geome- 
tries: 7"i = 5 /im, r 2 = 5.25 /mi. (B) a^ r — cr^ r , upper curves and the opposite a s rr — <r£ r , 
lower curves, calculated at the boundary zone between the semispherical apical and the 
cylindrical distal part of the growing pollen tube, parametrized by the wall thickness (up- 
per plots): r-i = 5.1 fim (green), r 2 = 5.2 /im (blue) and r 2 = 5.3 /tm (violet). The inner 
wall radius: r\ = 5 /im; turgor pressure P = 0.3 MPa. 
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Figure 5: Frustration energy Ep splitting due to topological effects (Insets: bifurcation 
diagram and tensegrity force F = — V£/(r) balance diagram) in a quasi-2D biological 
system - pollen tube apical region. Corresponding symmetry exchange takes place between 
the resonating residual energy levels E_ ± SE and E + ± SE of different major symmetry. 
Calculation performed at the transition zone between the (hemi-spherical) apical and the 
(cylindrical) distal part of a growing pollen tube (see Eq. (7)), at a constant turgor 
pressure P — 0.3 MPa. The inner and the outer wall radius in both subsystems read: 
ri = 5 /im, r2 = 5.25 /jm (wall thickness ~ 250 nm), respectively. The dispersion of each 
energy level SE = 0.0003. 
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Figure 6: Two branches of the anharmonic potential energy (frustration potential) 
U±(r) oc±4t4, possibly leading to pollen tubes growth rate oscillations, as a function 
of wall constituing molecules separation r. (Here: a = = 1; in general the coefficients 
a and /3 are linear in P, see Appendix 1). Oscillations take place between J7_ and U + 
potential energy level (by tunneling through symmetry change) yielding u), while the am- 
plitude is determined by the actual pressure P level, in accord with experiments (Kroeger 
and Geitmann, 2011a; Kroeger et al., 2011). 
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Figure 7: Least square fit of the experimental data (to) as a function of turgor pressure 
P for hypertonic (25 mM NaCl), isotonic and hypotonic (hypo-osmotic stress induced by 
the addition of water to the gel cultures) treatment of Nicotiana tobaccum pollen tube 
(Haduch and Pietruszka, 2012) fitted to the square root function (Appendix 1) derived in 
this paper (oj — 2irAy/P; [A 2 ] = m/kg). Stable turgor values correspond to those ranging 
between 0.1 and 0.4 MPa, which has been recorded using a turgor pressure probe, (Benkert 
et al., 1997). The initial point (0,0) added 'by hand' in the chart is exact - the pollen 
tube will not oscillate (w = 0) for the vanishing turgor pressure (P = 0). 
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